Comment on "Hole digging in ensembles of tunneling molecular magnets" 
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Tupitsyn et al. [Phys. Rev. B 69, 132406 (2004)] have recently reported results for the relaxation 
of crystalline systems of single-molecule magnets, such as Peg. They claim that, quite generally, 
(1) the magnetization and hole widths of field-distributions evolve with time t as \/i, and (2) the 
holes' line shapes are Lorentzian. We give a counter-example to these conclusions, and show that 
the main assumption on which some of them rest is invalid. 
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PIG. 1: The field function /(C,t) versus f/^^ for SC lattices 
and the shown values of the dipole field ^. The tunnel win- 
dow's half-width ^o is 0.1. The initial value of the magneti- 
zation 7710 is 0.2 of saturation. All data points are for systems 
of 32768 spins, and follow from averages over 800 runs, over 
times much greater than t ^ to but within the range where 
m oc t^. 

Tupitsyn, Stamp, and Prokof'eAii (TSP) have recently 
reported results for the relaxation of crystalline systems 
of single-molecule magnets, such as Feg, whose tunnel- 
ing windows are much narrower than their dipolar field 
spread. They find that, independently of crystalline lat- 
tice structure, (1) the magnetization m and hole widths 
of field-distributions evolve with time t as Vt, and (2) 
field-distribution line shapes are Lorentzian. The first 
conclusion is contrary to our prediction, of Ref. 2^, that, 
in zero field, m relaxes from weakly polarized states as 
tP, where p depends on crystal structure. We give below 
a counter-example to the conclusions TSP have reached, 
and show that the main assumption on which some of 
their conclusions rest is incorrect. 

We mainly use the notation of Ref. jj , giving the bias 
field S, in terms of the tunnel window field ^o i but ^o is 
given in terms of the nearest neighbor dipolar field Eo ■ 
We assume spins flip at rate I/tq if | ^ |< 1, but not 
at all otherwise, and time t is given in terms of tq. The 
following numbers may be found useful: the rms value of 
the dipolar field S^ is 3.7Ed and S,.3Ed for SC and FCC, 
respectively. 

Let P|(C,0 [Pi(^i'')] be the number density of up- 
spins (down-spins) with a field ^ acting on them, and 
let f{^,t) = piitt) - p^{tt). Note that m{t) = 




PIG. 2: (a) /(^, t) versus ^/t'' for PCC lattices and the shown 
values of t, ^o ~ 0.01, and p = 0.5. Lines are guides to the 
eye. (b) Same as in (a) but for ^o = 0.1. (c) Same as in 
(a) but for p — 0.726. The dashed line is for the best fitting 
Lorentz curve. The full line is from our theory, (d) Same as 
in (c) but for ^o ~ 0.1. In all cases, mo — 0.2 All data points 
are for systems of 65536 spins, and follow from averages over 
800 runs over times much greater than i 3> tq but within the 
range where m oc t^. All the shown data points come from the 
2^0 < C < 0.25^ (nearly one and two decades for ^o = 0.1 and 
^0 = 0.01, respectively). In the shown range of time values, 
1 — m/mo changes by approximately one and two decades (see 
Pig. 1), for 5o ~ 0.1 and ^o ~ 0.01, respectively. 



~ J d^f{(,,t). The main ingredient underlying Eq. 
(1) of Ref. [l| is the assumption that f{S,,t) ex 
7V(^) exp[— i/T(^)], where N{£^) is of no interest to us 
here, and t(^) is some time that depends only on S,- The 
Monte Carlo (MC) results shown in Fig. 1 are contrary 
to the assumption of TSP, that f{£,,t) is exponential in 
ti^ [The probability density that a spin have field £, and 
has not yet flipped at time "t" behaves much as /(^,t).] 

From the assumption that /(^, t) oc exp[— i/r(^)] and 
the further general statement TSP make, that l/r(^) is 
a Lorentzian function of ^, hole line widths that grow as 
^/t when t > To follow in Ref. 1] . Such ^/t growth does 
take place in SC lattices, but not in general, as one can 
gather from the MC results exhibited in Figs. 2a and 
2b for FCC lattices. Rather, from Figs. 2c and 2d, one 
gathers that ^ then scales as t^, where p ~ 0.73.^ Finally, 



in the region of interest, when t ^ tq but t is still within 
the range where m oc i^, the line shape in Figs. 2c and 
2d differ significantly from a Lorentzian function. For 
comparison, results from our own theorji^ are also shown 
in Figs. 2c and 2d. 
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FIG. 3: 1 — m/nio versus time t. All data points follow from 
MC simulations. Triangles are for SC lattices. All other data 
points are for FCC lattices. All data points for the FCC 
lattice are for systems of 8192 spins, except for •, which stand 
for 65536 spins. Data points for SC lattices are for 4096 spins. 
For all data points, averages over 1200 runs were performed, 
except for the 65536 spin system and t > 1200. For them, 
averages over 100 runs were performed. 

For completeness sake, we examine further numerical 
evidence that supports our claim that p = 0.73, not 1/2, 
for FCC lattices. To this end, note first that £, ^ t^ scal- 
ing implies, through the relation m{t) = — J d£_f{£_,t), 
that m ^ tP. Thus, the value of p can also be obtained 
from the time evolution of m. 

The difference between the relaxation of the magne- 
tization in SC and FCC lattices can be clearly appre- 
ciated in Fig. 3, as well as in Figs. 4a and 4b. Note 
in Fig. 3 that, for FCC lattices, m oc t^ and p ~ 0.73 
for time spans that are increasingly larger for larger val- 
ues of 1/^0- (This is as predicted in Ref. Q].) Note 
also that the p — 1/2 slope that is claimed by TSP to 
hold universally for all lattices appears to ensue for FCC 
lattices only when the relaxation crosses over from the 
m oz fP regime to saturation. This is in clear contrast 
with the data points for SC lattices in Fig. 3, for which 
p ~ 0.5. Linear m/rriQ versus time plots for SC and FCC 



lattices which further illustrate this point are also shown 
in Figs. 4a and 4b. The small neighborhood of the in- 
flection point (marked as a straight line segment) in Fig. 
4b is, of course, nearly straight. This transient behavior 
might be misinterpreted as the onset of a to ex \/t regime 
if, as in Fig. 1 of Ref. Q, data points below rajniQ ~ 0.5 
are not included in the plot. 

Finally, comparison of the results shown in Figs. 2a 
and 2b, Figs. 2c and 2d, as well as among different curves 
shown in Fig. 3 for different values of £o should allay any 
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FIG. 4: (a) ra/mo versus \/t. Data points are from averaging 
over at least 4000 MC runs for 4096 spins on SC lattices with 
^0 = 0.1. The straight line is a guide to the eye. (b) Same as 
in (a) but for 8192 spins on an FCC lattice. The straight line 
segment covers a neighborhood of the inflection point. 



misgivings about spurious effects that might arise from 
the finite number of spins Ug in the tunnel window, since 
Us ex fo if fo ^ ^£,- To the same end, data points for two 
FCC lattice sizes are shown in Fig. 3 for ^o — 0.01. 
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This is because df{^,t)/dt / —f{h,t)/r{^), since 
df{^,t)/dt depends, not only on / and ^, but also on the 
rate of magnetization reversal which is going on within the 
tunneling window at time t, which in turn depends non- 



trivially on time. else. 

This number is as given in our Phys. Rev. B paper in Ref. I. S. Tupitsyn, P. C. E . Stamp, and N. V. Prokof'ev, 

0. |arXiv:cond-mat/0407713l vl 27 Jul 2004. 
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